Enhancement of Dark Matter Annihilation via Breit-Wigner Resonance 
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The Breit-Wigner enhancement of the thermally averaged annihilation cross section (crw) is shown 
to provide a large boost factor when the dark matter annihilation process nears a narrow resonance. 
We explicitly demonstrate the evolution behavior of the Breit-Wigner enhanced {crv) as the function 
of universe temperature for both the physical and unphysical pole cases. It is found that both of 
the cases can lead a large enough boost factor to explain the recent PAMELA, ATIC, and PPB- 
BETS anomalies. We also calculate the coupling of the annihilation process, which is useful for an 
appropriate model building to give the desired dark matter relic density. 
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Introduction. The existence of dark matter is by now 
well confirmed [H. However, there is no candidate for the 
dark matter in the standard model. Understanding the 
nature of dark matter is one of the most challenging prob- 
lems in particle physics and cosmology. The recent cos- 
mological observations have established the concordance 
cosmological model where the present energy density con- 
sists of about 73% dark energy, 23% dark matter, and 4% 
atoms Currently, many dark matter search experi- 
ments are under way. These experiments can be classi- 
fied as the direct dark matter searches and the indirect 
dark matter searches. The direct dark matter detection 
experiments may observe the elastic scattering of dark 
matter particles with nuclei. The indirect dark matter 
searches are designed to detect the dark matter annihila- 
tion productions, which include neutrinos, gamma rays, 
electrons, positrons, protons and antiprotons. In addi- 
tion, the CERN LHC searches are complementary to the 
direct and indirect dark matter detection experiments. 

Recently, the indirect dark matter detection experi- 
ment PAMELA 3 reported an excess in the positron 
fraction from 10 to 100 GeV, but showed no excess for 
the antiproton data. The ATIC [l] and PPB-BETS [i| 
balloon experiments have also seen the excess in the 
6+ -I- e~ energy spectrum between 300 and 800 GcV. It is 
a natural idea that the dark matter annihilation can ac- 
count for the PAMELA, ATIC, and PPB-BETS anoma- 
lies. However, the thermally averaged annihilation cross 
section [uv) obtained from the observed relic density is 
far smaller than the required value from the PAMELA 
data. Therefore, one must resort to the large boost factor 
(about 100—1000) to explain the large positron flux. Cur- 
rent analysis on the dumpiness of dark matter structures 
indicates that the most probable boost factor should be 
less than 10 — 20 0. Considering the difficulty to yield a 
large boost factor, many authors investigate the decaying 
dark matter 0. However, the PAMELA data require the 
lifetime of dark matter to be of the order of 10^^ s. An 
alternative opinion is the nonperturbative Sommerfeld 
enhancement, which may provide a large boost factor as 
the weak force enhances the annihilation cross sections 



in the galactic halo [8|. 

The thermally averaged annihilation cross section {av) 
is a key quantity in the determination of the cosmic relic 
abundances of dark matter. On the other hand, {av) also 
determine the dark matter annihilation rate in the galac- 
tic halo. The only difference among the above two cases 
is the temperature T. For the relic density, {av) is usually 
evaluated at the freeze-out temperature x = m/T w 20 
(the averaged velocity v ~ yj'i/x), where m is the dark 
matter mass. The dark matter annihilation in the galac- 
tic halo occurs at a; « 3 x 10® {v « lO^'^). For nonrela- 
tivistic gases, {av) can usually be expanded in powers of 
X, {av) oc 0. For the s-wave annihilation {k = 0), 
{av) is a constant, which is independent of the tempera- 
ture of the Universe. For the p-wave annihilation [k = 1), 
{av) will be decreased as the Universe evolution. Clearly, 
only if fc < 0, {av) could be enhanced at the lower tem- 
perature. In such a case, one may obtain a large boost 
factor to explain the PAMELA, ATIC, and PPB-BETS 
anomalies. It is interesting to notice that when consider- 
ing the annihilation cross section at a narrow resonance, 
we can derive a negative number for k, which then indi- 
cates a Breit-Wigner enhancement mechanism. Recently, 
such an enhancement has explicitly been analyzed in Ref. 
[l3l (for the previous discussions, see Ref. [ll|). In the 
past, many authors have studied the dark matter anni- 
hilation near a resonance 

In this paper, we try to further give a comprehensive 
analysis on such a Breit-Wigner enhancement. Instead 
of using the center of mass frame, we work in the cos- 
mic comoving frame and adopt the usual single-integral 
formula to calculate {av). Except for checking the un- 
physical pole case, we will pay attention to the investiga- 
tion for the physical pole case in which the cross section 
{av) is found to have a maximum. In both cases, {av) 
will approach a constant as the Universe evolution. In 
terms of the observed dark matter abundance, we calcu- 
late the coupling of the annihilation process for the whole 
resonance parameter space. Hence, we derive the exact 
boost factor and find that both cases can lead a large 
enough boost factor to account for the PAMELA, ATIC, 
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and PPB-BETS results. 

Breit-Wigner enhancement. The PAMELA experi- 
ment observing no excess for the antiproton data indi- 
cates that dark matter wiU dominantly annihilate into 
the leptonic final states. In fact, the dark matter may 
first annihilate into some particles, such as the Higgs 
triplets in the left-right symmetric model [l^l, and then 
these particles decay into the charged leptons. For the 
purpose of this paper, we simply consider that two dark 
matter particles directly annihilate into a pair of charged 
leptons via S channel Higgs boson exchanging. Since the 
Breit-Wigner enhancement requires a narrow resonance, 
we follow Ref. [lo| to introduce an auxiliary parameter 
S (I (5 1 ^ 1) to express the intermediate particle mass M 



Ar 



4m^(l - S) 
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For the S < case, one may obtain a physical pole. In 
this case, the exchanging particle may decay into both 
initial states and final states. For a given decay width F, 
one may write the following annihilation cross section 
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where Bi and Bf are the branching fractions of the 
resonance into the initial and final channels, respec- 
tively. Because oi Bi + B f = 1, one can directly obtain 
BiBf < 0.25. Here we have neglected the masses of final 
leptons. The parameter s is defined by s = {pi -f ^2)^, 
where pi and p2 are the four-momenta of initial dark 
matter particles. For the S > case, we have an unphys- 
ical pole. In this case, the intermediate particle can not 
decay into the initial dark matter particles. Therefore, 
we introduce a vertex a for the trilinear coupling among 
two dark matter particles and the exchanging particle to 
express the annihilation cross section 



4:EiE2<7v = 2a^ 
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For the thermally averaged annihilation cross section 
(av) , we adopt the usual single- integral formula 




FIG. 1: Numerical illustration of the Breit-Wigner enhanced 
(av) as a function of x. The parameter Ro is defined as 
Ro = {av) / {(tv)q, with {av)o denoting the thermally aver- 
aged annihilation cross section at T = 0. 
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FIG. 2: Numerical illustration of the Breit-Wigner enhanced 
(cru) as a function of x. The parameter Rj is defined as 
Rf = {crv) / {av)x^20, with {av):s=2Q denoting the thermally 
averaged annihilation cross section at the usual "freeze-out" 
time X = 20. 
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where Ki{x) and K2{x) are the modified Bessel func- 
tions, gi = 1 is the internal degrees of freedom of dark 
matter particle. Since v of {av) is the M0ller velocity, 
we work in the cosmic comoving frame. If one takes the 
center of mass frame [loj , {av) should be multiplied by a 
factor {1 + Kl{x)/Kl{x))/2 



For an analytical illustration, we redefine 

s = 4to2(1 + z). (7) 

Then the integration region of Eq. ^ becomes < z < 
00. The annihilation cross sections AEiE2av in Eqs. ([2]) 
and (BI) can be rewritten as 



AEiE2av (X 



l + z 
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where 7 is given by 



7 = F/M 



(8) 
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For 2; > 20 and \S\,j < 0.1, {av) in Eq. ^ can approxi- 
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mately be given by 



(trw) (X x'' 



dz 
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It is worthwhile to stress that the integration result of Eq. 
pH)) is insensitive to x when z is negligible in (2+^)^+7^. 
In the 5 > Q case, the effective integration upper bound 
is Zeff ^ 4/a;, thus {av) can be enhanced as the uni- 
verse evolution. For the 5 < Q case, one may derive 
Ze// ~ niax[4/x, 2|(5|] when |^| > 7, we then find that 
{av) has a maximum at x ~ 2/|(5|. If |i5| ^ 7, one can- 
not obtain an obvious peak. When x 3> 4/max[|5|, 7], 
((Tw) will approach to a constant for both cases. Our 
numerical results [using Eq. (|4])] in Figs. [1] and [2] ex- 
plicitly demonstrate the above analysis. In Fig. [l] the 
parameter i?o is defined as Rq = {av) / {av)a, with ((tw)o 
denoting the thermally averaged annihilation cross sec- 
tion at r = 0. One may easily see from Fig. [T]that the 
i5 < case gives the larger Rq at the higher temperature. 
In Fig. [21 we plot the ratio Rj = {av) / {av)x=20- It is 
seen that for small \8\ and 7 with |5| ~ 7 < 10"'^, both 
of the cases provide a significant enhancement for the 
thermally averaged annihilation cross section. For much 
smaller = 7 = 10~^, the 5 < case can give the larger 
enhancement. It should be mentioned that our results 
are independent of a, BiBf, and m. 

Boost factor and couplings. The Breit-Wigner en- 
hancement mechanism can change the indirect dark mat- 
ter detection. On the other hand, the Breit-Wigner en- 
hancement will affect the calculation of the dark matter 
relic density because the annihilation process does not 
freeze out even after the usual "freeze-out" time Xf — 20 
p^ . More importantly, we can derive the size of a 
{6 > 0) and BiBf {S < 0) for given 5 and 7 from the ob- 
served dark matter abundance Vtuh? = 0.1131 ± 0.0034 
0. The values of a and BiBf can help us to build ap- 
propriate models. With the help of Eqs. ([2]), ([3]) and 
(ID), one can calculate the thermally averaged annihila- 
tion cross section {av) in the galactic halo, which allows 
us to obtain the exact boost factor. It is worthwhile to 
stress that the Breit-Wigner enhancement docs not affect 
the direct dark matter searches. 

The evolution of dark matter abundance is given by 
the following Boltzmann equation Q: 



dY 
dx 
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where Y = n/s{x) denotes the dark matter number den- 
sity. The entropy density s(a;) and the Hubble parameter 
H evaluated at a; = 1 are given by 



s(x) 
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45 MpL 



(12) 
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where Mp^ — 1.22 x 10^^ GeV is the Planck energy, 
is the total number of effectively relativistic degrees of 
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FIG. 3: Numerical illustration of the boost factor BF (solid 
lines) and the coupling a (dashed lines) on the S and 7 planes 
for the S > case. 
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FIG. 4: Numerical illustration of the boost factor BF (solid 
lines) and the parameter BiBf (dashed lines) on the 5 and 7 
planes for the 5 < case. 



freedom. Here we choose = 106.75 for illustration. 
Using the result Yq of the integration of Eq. pT|) . we 
may obtain the dark matter relic density f^^i/i^ 



nnh^ = 2.74 X 10*^ 



GeV 



:Yn 



(14) 



Using the Boltzmann equation in Eq. pT|) . we nu- 
merically calculate a and BiBf for the unphysical pole 
case and the physical pole case, respectively. Our nu- 
merical results (dashed lines) are shown in Figs. [3] and 
m Here we have taken to = 1 TeV. For the S > case, 
we obtain 15 GeV ^ a < 4.9 TeV, which can be easily 
satisfied. For the S < case, the parameter BiBf is far 
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less than the upper bound 0.25 except for the lower left 
region. For most of the parameter range, B^Bf < 0.001 
indicates that the successful models must have a hierar- 
chy between the initial branching factor Bi and the final 
branching factor Bf. If the dark matter mass m is en- 
larged by N times, a and BiBf in Figs. [3] and |4] should 
be approximately enlarged by iV^ times. 

After obtaining a and BiBf, we calculate the ther- 
mally averaged annihilation cross section (av) in the 
galactic halo (x « 3 x 10^), which is shown in Figs. [3]and 
131 We would like to emphasize that our results are in- 
sensitive to the dark matter mass m. Here we normalize 
(av) by the usual nonresonance annihilation cross section 
10"^ GeV"^ to define the boost factor BF 



BF 



10-9 GeV" 



(15) 



It is clear that smaller \6\ and 7 will provide larger boost 
factors. For the S > case, the large boost factor {BF > 
100) requires 5, 7 < 0(10"^). Our results have some 
differences from Fig. 4 in Ref. [l3| for the region BF < 
C(10) even if we choose = 200. For the S < case, one 
may obtain (5,7 < 0(10"'') for BF > 100. In the lower 
left region of Fig. |4l we find BF <^ 1, which implies 
that the indirect dark matter detection experiments will 
not find any signal of the dark matter annihilation. It is 
clear that both cases can provide a large enough boost 
factor to explain the PAMELA, ATIC, and PPB-BETS 
anomalies. This is one of our primary results. 

Discussion and Conclusion. We have evaluated the 
boost factor BF and the values of a, BiBf in terms of the 
observed dark matter abundance fljjh^. If the parame- 
ter or BiBf is enlarged by N times, the Breit-Wigner 
enhanced annihilation cross section will be enlarged by 
the same times. While the dark matter relic number 



density Y will be approximately suppressed by N times, 
one thus needs to introduce new dark matter candidates. 
Although we can obtain the larger boost factor, this sce- 
nario will give the smaller dark matter annihilation rate, 
which is proportional to (av) Y^. In fact, many mod- 
els have several dominant annihilation processes, which 
may include the nonresonance and resonance cases. If the 
nonresonance annihilation processes determine the dark 
matter relic density, one cannot obtain a, BiBf, and BF 
from fliyh^. In this case, one may determine those pa- 
rameters from other constraints; the required values for 
a, BiBf, and BF must be smaller than the predicted 
values in Figs. [3]and[4l In general, it still has some pa- 
rameter space in this case to account for the PAMELA, 
ATIC, and PPB-BETS resufts. 

In conclusion, we have made a comprehensive analysis 
based on the Breit-Wigner enhancement near the reso- 
nance point. In terms of the observed value oi ilnh^, we 
have evaluated the exact thermally averaged annihilation 
cross section {av) in the galactic halo and the boost factor 
BF, and calculated the couplings a and BiBf of the an- 
nihilation process for both the 6 > and the 6 < cases, 
respectively. The numerical results lead us to a general 
conclusion that both the (5 > and the 6 < cases can 
provide a large enough boost factor BF > 100 to ex- 
plain the PAMELA, ATIC, and PPB-BETS anomalies. 
It would be interesting to find a model with the appro- 
priate coupling a or BiBf for the annihilation processes 
to give the desired dark matter relic density. 
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